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O ; Abstract 

Ph. We construct a simple 5*0(10) SUSY GUT with D 3 family symmetry and 

Q-t, low energy R parity The model describes fermion mass matrices with 14 

parameters and gives excellent fits to 20 observable masses and mixing an- 
gles in both quark and lepton sectors, giving 6 predictions. Bi-large neutrino 
mixing is obtained with hierarchical quark and lepton Yukawa matrices; thus 
avoiding the possibility of large lepton flavor violation. The model naturally 
predicts small 1-3 neutrino mixing, sin #13 ~ 0.05, and a CP violating phase 
S close to 7r/2. Among other interesting predictions is a tiny effective Majo- 
rana mass for neutrinoless double-beta decay. Leptogenesis is also possible 
with the decay of the lightest right-handed neutrino giving an acceptable CP 
violating asymmetry e x of order 10~ 6 , and with the correct sign for the resul- 
tant baryon asymmetry. Note, similar models with the non-abelian symmetry 
groups SU(2) or _D 4 , instead of D 3 , can be constructed. 

There are very few experimental indications for new physics beyond the extremely 
successful standard model. By far the two most exciting harbingers of new physics are 
the successful prediction of gauge coupling unification in supersymmetric GUTs and 
the discovery of neutrino masses and mixing. Both of these experimental results hint at a 
new large mass scale of order 10 16 GeV. In addition, SUSY GUTs can alleviate the gauge 
hierarchy problem implying the exciting possibility for the discovery of supersymmetric 
partners at the LHC. It is thus important to see to what extent SUSY GUTs and neutrino 
masses and mixing angles are compatible. 



The first problem one encounters in any theory with quark-lepton unification is the 
fact that flavor mixing in the quark sector is small, while for neutrinos there are two large 
mixing angles. In the quark sector it is well-known that hierarchical Yukawa matrices 
"naturally" fit the family hierarchy, as well as the hierarchy of flavor mixing evident in 
the CKM matrix. Moreover this can be described group theoretically in terms of the 
hierarchical breaking of family symmetries via the Froggatt-Nielsen mechanism [2] . In 
super symmetric theories, family symmetries play a dual role. In addition to describing 
the hierarchy of fermion masses and mixing, they "naturally" align squark mass matrices 
with their quark superpartners; 1 thus ameliorating problems with flavor violation. In 
this letter, we show that bi-large neutrino mixing is naturally incorporated into SUSY 
GUTs with hierarchical quark and lepton Yukawa matrices; thus avoiding the possibility 
of large lepton flavor violation 2 . Given Yukawa matrices of the form 




Y~\baa\\ (1) 
\ b a 1 / 

with b <C a <C 1 and a hierarchical right-handed neutrino mass matrix of the form 

/ M Rl \ 
M R =\ M R2 (2) 
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with M Rl <C M R2 -C M Rs one "naturally" obtains large neutrino mixing [7J |Hl EH EH EJ 
H2HH]. Neglecting, for example, the contribution of M R23 we find, upon integrating out 
the lightest right-handed neutrino 

m.= \ H M^J! (3 ) 

giving maximal mixing angle for atmospheric neutrinos. It has been shown by one of us 
that this mechanism also gives acceptable bi-large neutrino mixing if the contributions 
of the lightest two neutrinos are comparable [Tlj . 

Our model is an SO(10) realization of this general mechanism. This letter is organized 
as follows: first we present the superpotential for the fermion mass sector of an 5*0(10) 

1 For non-abelian family symmetries, squark and slepton mass matrices are proportional to the unit 
matrix prior to family symmetry breaking. Then when the family symmetry is broken, generating the 
hierarchy of fermion masses, the squark and slepton mass matrices also receive corrections which are 
aligned with their fermionic partners. 

2 In several recent SUSY GUT constructions, bi-large neutrino mixing is obtained with so-called 
lopsided Yukawa matrices for quarks and leptons [HJ- In this approach, use is made of the SU (5) relation 
Y e = Yj where large, unobservable right-handed mixing for down quarks is converted to large left- 
handed mixing angles for leptons (for a review, see 0JE]). In this case, one takes (Yd) 23 ~ (Yd) 33 ~ 1- 
One apparent difficulty of this approach is the large induced lepton flavor violation for processes such 
as /i — > e7 due to large off-diagonal (Y e ) 32 « (Y e ) 33 ~ 1 0. 



2 



SUSY GUT. We then fit the low energy data using a global x 2 analysis. The y 2 function 
includes 20 low energy observables for fermion masses and mixing angles, in addition to 
the three low energy gauge couplings. We obtain an excellent fit to the data. Then, we 
present predictions for additional observables (not included in the x 2 analysis): sin 2 #13, 
the CP violating angle 5 (and Jarlskog parameter J), the effective neutrinoless double 
beta decay mass parameter (m)^, the tritium beta decay mass parameter m e JJ ', and 
finally the lepton number asymmetry parameter relevant for leptogenesis, e±. Alternate 
versions of the model with family symmetries SU (2) or _D 4 (instead of D 3 ) can easily be 
constructed [T3] . 

The model 

Consider an 50 (10) SUSY GUT with a D 3 x [U(l) x Z 2 x Z 3 ] family symmetry. 3 
The three families of quarks and leptons are contained in three 16 dimensional repre- 
sentations of 50(10) {16 a , I63} with 16 a , a — 1,2 a D3 flavor doublet (see appendices 
of Ref. [T^] for details on D3). Consider the charged fermion sector first. Although 
the charged fermion sector is not the main focus of this letter it is, however, necessary 
to first construct a superpotential which is consistent with charged fermion masses and 
mixing angles, since in 50(10) the neutrino and charged fermion sectors are inextricably 
intertwined. 

The superpotential resulting in charged fermion masses and mixing angles is given 

by 

.f ermions I63 10 I63 + 16 a 10 X a (4) 

+Xa (M x X a + 45 f| 16 3 + 45 § 16„ + A 16„) 

where M x = Mq(1 + aX + f3Y) includes 5O(10) breaking vevs in the X and Y directions, 
a , (fi a (D 3 doublets), A (1b singlet) are 5O(10) singlet flavon fields, and M, M are 
5O(10) singlet masses. The fields 45, A, 0, are assumed to obtain vevs (45) ~ 
(B - L) M G , A < M and 

<*>-(£)•<»-(£) (5) 

with 0i > 02. 

The superpotential, (Eqn. Hj) results in the following charged fermion Yukawa matri- 

3 For related charged fermion analyses in 5O(10) SUSY GUTS with D 3 x U{1) (or SU{2) x U{l) n ) 
family symmetries, see ^H] (or ^fOElE])- Note, the novel feature in this letter is the structure and 
results for neutrino masses and mixing. 
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Y d = -e' i -ea A (6) 



Y e = e' 3 e 3 e | A 



with 

£ = 2 /0 i; e ex 2 /M; (7) 
e oc &/M; e' ~ (A/M ); 
1 + a 

a = — ; p ~ < a. 

1 — da 

Let us now discuss neutrino masses. In the three 16s we have three electroweak 
doublet neutrinos (u a , v 3 ) and three electroweak singlet anti-neutrinos (u a , P 3 ). 5 The 
superpotential W c h,f ermicms also results in a neutrino Yukawa matrix: 




Y v = ( e' u 3 e w § e | A (8) 

with = 2a /{2a — 1) and a Dirac neutrino mass matrix given by 

v 

m u = Y u —= sin p. (9) 
V2 

In addition, the anti- neutrinos get GUT scale masses by mixing with three 5*0(10) 
singlets {N a , a = 1,2; N 3 } transforming as a D 3 doublet and singlet respectively. The 
full superpotential is given by W = W c h .fermicms H~ ^^neutrino with 

W neutrino = 16 (A 2 N a 16 a + A 3 iV 3 16 3 ) (10) 
+ | (5 a iV a N a + S 3 N 3 N 3 ) . 

We assume 16 obtains a vev, w 16 , in the right-handed neutrino direction, and (S a ) = M a 
for a = 1, 2 (with M 2 > M x ) and (S 3 ) = M 3 . 6 

4 In our notation, Yukawa matrices couple electroweak doublets on the left to singlets on the right. It 
has been shown in Ref. |17j that excellent fits to charged fcrmion masses and mixing angles are obtained 
with this Yukawa structure. 

In an equivalent notation, we have three left-handed neutrinos (vta = v ai v ls = ^3) and three 
right-handed neutrinos defined by (vR a = ^, vrs, = ^3)- 

6 These are the most general set of vevs for <p a and S a ■ The zero vev for (pi can be enforced with a 
simple superpotential term such as S 4> a 4>a- 
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We thus obtain the effective neutrino mass terms given by 



W = vm u v + vVN+-NM N N (11) 

2 

with 

/ A 2 \ 

V = v 16 A 2 , Mjv = diag{M x , M 2 , M 3 ). (12) 
V A 3 / 

A simple family symmetry giving the desired form of the superpotential 7 is D 3 x 
U(l) x Z 2 x Z 3 where the -D 3 charges were defined earlier, while the U(l) charge as- 
signments are (1 for 16 3 , 2 for 16 a , -2 for N a , -1 for A^, -1 for 45, for 16 and x a ) an d 
everyone else fixed by these. In addition we assign Z 2 charges (16 3 , 16 a , N 3 , N a , x a , Xa) 
odd, all others even and Z 3 charges a = e"^ for all fields, except 45 with charge 1. 
Note, that Z 2 can also be interpreted as a family reflection symmetry which guarantees 
an unbroken low energy R parity [TH] . 

The electroweak singlet neutrinos {p, N} have large masses V, Mn ~ Mq. After 
integrating out these heavy neutrinos, we obtain the light neutrino mass matrix given 
by 

M =m v M^ 1 ml, (13) 
where the effective right-handed neutrino Majorana mass matrix is given by: 

M R = V M^ 1 V T = diag(M fll , Mr 2 , Mr 3 ), (14) 

with 

M Rl = (\ 2 v 16 ) 2 /M 2 , Mr 2 = {\ 2 v 16 ) 2 /M 1 , Mr 3 = (\ 3 v 16 ) 2 /M 3 . (15) 

Defining U e as the 3x3 unitary matrix for left-handed leptons needed to diagonal- 
ize Y e (Eqn. EJ, i.e. Y e D = Y e U* E and also U v such that C/J M U v = M D = 
diag(m vi , m„ 2 , m V3 ), then the neutrino mixing matrix is given by Upmns — U\ U v in 
terms of the flavor eigenstate (u a , a = e, fi, r) and mass eigenstate (z/j, i — 1,2, 3) basis 
fields with 



^2(UpMNs)ai Vi- (16) 



For Upmns we use the notation of Ref ED] with 









H 







ci 2 ci 3 s 12 c 13 s 13 e iS \ / e iai / 2 vi 

-si 2 c 23 - c 12 s 23 si 3 e i<5 ci 2 c 23 - si 2 s 23 si 3 e i<5 s 23 ci 3 e ia2/2 u 2 | (17) 

iS iS I \ 

Si2 s 23 — Ci 2 c 23 Si 3 e — Ci 2 s 23 — Si 2 c 23 Si 3 e c 23 Ci 3 / \ v 3 



7 Note, highly suppressed terms of the form S a N a N3 (j> a §\ are still allowed. This too can be 
forbidden by an additional discrete Z4 symmetry. 
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The 3x3 Majorana mass matrix is of the general form discussed by many authors 
j. We have 

M =Vi Mi Vi +V 2 M 2 V2 + V3 M 3 V 3 (18) 

with 

/ \ /\ ■ r\ 2 

Mi= (e'M) 2 -3e'e\iau\ ( A " ) M 2 ; (19) 

\0 -3e' e\iau\ (3 e |f a\f ) W2 A 2 W 

jW 2 = -3e'\~euj 2 \ (3|eo;|) 2 -9 e |e a uo\ P Mi; (20) 

VSe' e|a^| -9e|eacu| (3 e |a|) 2 / ^2 A 2 v 16 J 



-3 
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and 



M»=| lel(|e|^|) 2 (| ^ |c^|) 2 fcH I ( "\ ) M 3 (21) 



p x = e^*" ; (22) 
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p 2 = e -*(*e+*»J ; (23) 




V 3 = e^*" . (24) 



■Mi, -A/f 2 , A^3 are in general complex rank 1 mass matrices. However only the difference 
in their overall phases may be observable. Thus, there are, in principle, two new CP 
violating phases in the neutrino sector, in addition to the four phases already fixed by 
charged fermion masses and mixing angles. We shall impose the constraint that neutrino 
Majorana masses Mi are all real. This eliminates two arbitrary phases. We note that 
the best fits with free phases for Mj are very close numerically to our zero phase model. 

Fitting the low energy data : x 2 analysis 

Yukawa matrices in this model are described by seven real parameters {A, e, e, a, 
p, e', £} and, in general, four phases $ p , $5}. Therefore, in the charged 

fermion sector we have 11 parameters to explain 9 masses and three mixing angles 
and one CP violating phase in the CKM matrix, leaving us with 2 predictions. 8 Note, 
these parameters also determine the neutrino Yukawa matrix. Finally, our minimal 

8 Of course, in any supersymmetric theory there is one additional parameter in the fermion mass 
matrices, i.e. tan [3. Including this parameter, there is one less prediction for fermion masses, but then 
(once SUSY is discovered) we have one more prediction. This is why we have not included it explicitly 
in the preceding discussion. 
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ansatz for the right-handed neutrino mass matrix is given in terms of three additional 
real parameters 9 , i.e. the three right-handed neutrino masses. At this point the three 
light neutrino masses and the neutrino mixing matrix, Upmns, (3+4 observables) are 
completely specified. Altogether, the model describes 20 observables in the quark and 
lepton sectors with 14 parameters, effectively having 6 predictions. 10 

In addition to the parameters describing the fermion mass matrices, we have to input 
three parameters specifying the three gauge couplings: the GUT scale M G defined as 
the scale at which a\ and a 2 unify, the gauge coupling at the GUT scale Qg, and the 
correction e 3 to a^Mc) necessary to fit the low energy value of the strong coupling 
constant. Finally we have to input the complete set of soft SUSY breaking parameters 
and the value of fi(Mz)- 11 All the parameters (except for /j,(Mz)) are run from the GUT 
scale to the weak scale (Mz) using two (one) loop RGEs for dimensionless (dimensionful) 
parameters. At the weak scale, the SUSY partners are integrated out leaving the two 
Higgs doublet model as an effective theory. We require proper electroweak symmetry 
breaking. Moreover, the full set of one loop, electroweak and SUSY, threshold corrections 
to fermion mass matrices are calculated at Mz- Below Mz we use 3 loop QCD and 1 loop 
QED RG equations to calculate light fermion masses. More details about the analysis 
can be found in or [21]. 12 In addition, we self-consistently include the contributions 
of the right-handed neutrinos to the RG running between the GUT scale and the mass 
of the heaviest right-handed neutrino |25j . 

The x 2 function is constructed from observables given in Table H Note that we over 
constrain the quark sector. This is due to the fact that quark masses are not known 
with high accuracy and different combinations of quark masses usually have independent 
experimental and theoretical uncertainties. Thus we include three observables for the 
charm and bottom quark masses: the MS running masses (m c (m c ), mb{mi,)) and the 
difference in pole masses Mf, — M c obtained from heavy quark effective theory. The same 
is true for observables in the CKM matrix. For example, we include Vu and the two CP 
violating observables €k and the value for sin(2/3) measured via the process B —>■ J/ip K s . 
We thus have 16 observables in the quark and charged lepton sectors. We use the central 
experimental values and one sigma error bars from the particle data group j20] • In case 
the experimental error is less than 0.1% we use a = 0.1% due to the numerical precision 
of our calculation. 

At present only four observables in the neutrino sector are measured. These are the 
two neutrino mass squared differences, Am^ and Am^, and two mixing angles, sin 2 # 12 
and sin 2 #23 ■ For these observables we use the central values and 2a errors from Ref. [26J. 
The other observables: neutrino masses, 1-3 mixing angle and the phase of the lepton 

9 In principle, these parameters can be complex. We will nevertheless assume that they are real; 
hence there are no additional CP violating phases in the neutrino sector. 

10 Note, the two Majorana phases are in principle observable, for example, in neutrinoless double-beta 
decay [21], however, the measurement would be very difficult (perhaps too difficult [22])- If observable, 
this would increase the number of predictions to 8. 

11 The reader will find more details on the soft SUSY breaking parameters in the next section. 

12 The only difference is that in the present analysis wc include all three families of fcrmions. 
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mixing matrix are predictions of the model. We will discuss these later. 
Best fit 

In Table [2 we present the best fit for quark and lepton masses and mixing for the \ 2 
function constructed from the 20 observables in the quark and lepton sectors and 3 gauge 
couplings. The input parameters for this fit are summarized at the top of the table. One 
can see that the charged fermion sector is fit very well with the largest contributions to 
X 2 coming from m^/m s and sin (2 ft). 13 

The neutrino observables, on the other hand, are fit very close to their central values, 
giving a negligible contribution to x 2 - Note, our model does not rely on single right- 
handed neutrino dominance. In the case of SRHND [3 El El EH EH E2j one obtains 
maximal atmospheric mixing angle as described (see Eqn. OJ), but these models (and 
others) differ significantly in the way they obtain the large solar mixing angle. Many 
models actually abandon the hierarchical texture of the neutrino Yukawa matrix or they 
assume a different form for the right-handed neutrino mass matrix. In [T^j it was 
shown that it is impossible to get the desired "mild" hierarchy in m V2 /m V3 when one 
of the right-handed neutrinos dominates in strictly hierarchical models. To fix this, for 
example in j2Zj ; the authors make the (1 2) element comparable to the (2 2) element. 
But the Yukawa matrix is no longer hierarchical as in Eqn. ^ It appears that the only 
possible way to get the correct neutrino masses, and at the same time have bi-large 
mixing with Yukawa matrices of the form in Eqn. ^ is for the contributions of Mr x 
and Mr 2 to be comparable (but with opposite sign). Note that the (2 3) block of the 
resulting M v does not change when varying the relative contribution of Mr x and Mr 2 . In 
[T%] . a good fit was found with Mr 1 contributing only -1.5 times more than Mr 2 . In our 
model, for the best fit, we find that Mr x contributes only -1.2 times more than Mr 2 . u 
In summary, comparable contributions of and Mr 2 are essential for getting bi-large 
mixing and an "observed" neutrino mass hierarchy (neutrino mass squared differences 
interpreted as normal mass hierarchy) in hierachical models with Yukawa matrices of the 
form Eqn. ^ Finally the contribution of Mr s is, in a model independent analysis, not 
very constrained, since we do not know the mass of the lightest neutrino. However, for 
the best fit in our model it turns out that all three right-handed neutrinos contribute 
comparably, and there is no single right-handed neutrino dominance. Indeed, forcing the 
contribution of M^ 3 to be negligible (by pushing it all the way to the GUT scale) makes 
the fit worse, but does not change the fact that we would have bi-large mixing. 

The best fit to fermion masses prefers the region of SUSY parameter space charac- 
terized by fx, Mi/2 *C mi6 an d — A ~ v^mio — 2m 16 , where M 1 / 2 is the universal 
gaugino mass, mio is the universal Higgs mass, mi 6 is the universal squark and slepton 

13 The problem associated with md/m s is directly connected with the Georgi-Jarlskog factor of 
three explicit in our model, while the problem with sin(2/3) is a result of using factorizable phases in 
our Yukawa matrices [2SI ■ S.R. thanks G.C. Branco for pointing out the origin of the latter problem. 

14 We define n = (\Y n {l,i)\ 2 + \Y n (2,i)\ 2 + \Y n (3, i)\ 2 )/M Ri as the contribution of Mr, to M v . This 
definition is not unique, there is no reason to take i)\ 2 , ... because 1^(1, i) 2 , ... is what appears in 

M v . But the latter is complex, so for a rough order of merit we take the absolute values instead. 
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mass and Aq is the universal trilinear coupling. 15 This is required in order to fit the top, 
bottom and tau masses when the third generation Yukawa couplings unify [22] • Note, 
third generation Yukawa unification receives only small corrections from 2-3 mixing in 
our model. For completeness, predictions for squark, slepton and Higgs masses for the 
best fit are presented in Table El However, it should be emphasized that the x 2 function 
is not very sensitive to changes in the SUSY parameters, as long as the relations discussed 
above are approximately satisfied. 16 As a consequence of these relations we expect heavy 
first and second generation squarks and sleptons, while the third generation scalars are 
significantly lighter (with the stop generically the lightest). In addition, charginos and 
neutralinos are typically the lightest superpartners. We predict values of tan /3 ~ 50 and 
a light Higgs mass of order 120 GeV. The specific relations between the SUSY break- 
ing parameters also leads to an interesting prediction for the process B s — > fi + fi~ with 
branching ratio in the region currently being explored at the Tevatron. 17 Furthermore, 
the neutralino relic density obtained for our best fit parameters is consistent with WMAP 
data j21] and direct neutralino detection is possible in near future experiments. Finally, 
this region maximally suppresses the dimension five contribution to proton decay [3U] 
and suppresses SUSY flavor and CP violation in general. For a detailed analysis of SUSY 
and Higgs spectra and related phenomenology see Ref. and [2~lj . 

Additional predictions 

Finally let us discuss the predictions in the lepton sector, summarized in Table El At 
this point we would like to emphasize, that non of these observables were included in the 
X 2 analysis and no parameters were constrained or tweaked in order to get the values 
we present. Therefore these observables are genuine predictions of the model. However, 
it should be noted that slightly different values of some of these observables might be 
obtained without significant changes in x 2 ■ We did not study in detail to what extent 
each of the observables can be modified without significantly worsening the fit. 

The only mixing angle in the lepton sector, which so far has not been measured, is 
1-3 mixing. The 2a upper bound from the global fit to neutrino data (with Chooz and 
solar+KamLAND data contributing comparably to the fit) is |26j : 



15 In addition, we require non-universal Higgs masses at the GUT scale. For the best fit we need 
l/2{m 2 Hd — r rn 2 H J)/m\ ~ .07. Note this is significantly smaller than needed in the past [231 • That is 
because the RGE running of neutrino Yukawas from Mg to the heaviest right handed neutrino has 
been included self-consistently. As noted in |2.'{| , such running was a possible source for Higgs splitting. 
Evidently it can not be the only source. 

16 Note, the three input parameters (/z, Mx/2, itiiq) are not varied when minimizing \ 2 . Moreover, 
X 2 is basically flat with increasing m\Q. 

17 This process is sensitive to the CP odd Higgs mass, rriA, which can be adjusted in theories with 
non-universal Higgs masses. 



sin 2 01, < 0.031. 



(25) 



We obtain the best fit value given by 




(26) 
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This is well below the upper bound and unfortunately too small to be seen in current or 
near future experiments. In particular, after three years of running the Double Chooz 
experiment will only probe down to sin 2 2^13 > 0.03 [2^ an d the maximal sensitivity 
goal for other future reactor experiments is just at the border of observability with 
sin 2 26*13 > 0.01 [S2]- On the other hand, long baseline accelerator experiments, such as 
NUMI and T2K, may be able to probe the best fit value |3*3*] . 

CP violation in the lepton sector given by, for example, the Jarlskog parameter, 

J = ImiUeiU^U^U/a) = s 12 c 12 s 13 c 2 13 S23C23 sin(<f), (27) 

(where we use the standard parametrization, Eqn. El with the abbreviated notation 
Sjj = sin^jj and c$j = cos 6V,), is proportional to sin^. In spite of the fact that the best 
fit suggests a large CP violating phase 5 ~ ir/2, the smallness of sin 6*13 results in J ~ 
0.01. CP violation of this magnitude may be observable at long baseline experiments. For 
example, the JPARC-SK experiment has a potential sensitivity to sin 2 2#i3 < 1.5 x 10 -3 
and 5 ~ ±20° and a comparable sensitivity is expected from the "Off-axis NUMI" 
proposal • 

Majorana neutrino masses are, in principle, observable in processes like neutrinoless 
double-beta decay. The effective mass 



2 



2 m„, e l < 



(28) 



(where a- = a« + 25, % = 1, 2 [35]) is predicted to be of order 2 x 10~ 4 eV which is too 
low to see in near- future experiments (UHl EH]. This is a consequence of the Majorana 
phases a± and a 2 being almost opposite (see Table EJ). The contribution of m Ul and 
m U2 is larger than m Vz by an order of magnitude and, due to the opposite Majorana 
phases, these tend to cancel. This is also evident from the prediction for the effective 
electron-neutrino mass observable, relevant for the analysis of the low energy beta decay 
of tritium. This mass parameter is unaffected by Majorana phases and is predicted to 
be an order of magnitude larger. The observable, 



1/2 



El^| 2 < ( 29 ) 



is predicted to be of order 6 x 10~ 3 eV. The current experimental limit is m e JJ < 2.5 
eV with the possibility of future experiments, such as KATRIN, reaching bounds on the 
order of 0.35 eV [33]. Unfortunately, both mass parameters may be unobservable by 
presently proposed experiments. 

The best fit values of the heavy right-handed neutrino masses, defined in Eqn. E3 
are given by 

M Rl , M R2 , Mr 3 w 10 10 , 10 12 , 7xl0 14 GeV . (30) 
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The lightest Majorana neutrino, Rx, is responsible for leptogenesis. The lepton number 
asymmetry parameter is given by |3T1 13*%] 



* =-Tfc(^feE^ 3 M([^]H) 2 ^} (31) 

where Y v is the Dirac neutrino Yukawa matrix 18 . This formula is only valid in the 
limit Mrj << Mr 2 ,Mr 3 . An acceptable baryon asymmetry, obtained via leptogenesis, 
requires values of e± ~ O(10~ 6 ). In comparison, for the relevant parameter for non- 
thermal leptogenesis we find 4 x e\ ~ —4.6 x 10~ 7 . Note that the baryon asymmetry 
obtained from leptogenesis comes via electroweak baryon and lepton number violating 
interactions preserving B - L. Hence the resultant baryon asymmetry Nb oc —Nl oc —t\. 
Hence we even obtain the correct sign. 19 

For the case of thermal leptogenesis, the relevant asymmetry parameter is 8 x e\ ~ 
—.9 x 10 -6 . 20 We note that the values for the parameters M Rl ~ 10 10 GeV and m Ul pa 
4 x 10~ 3 eV are in the regime where the cosmological baryon asymmetry only depends on 
neutrino properties j2Hj- 21 It does not depend on the initial neutrino abundance nor the 
initial baryon asymmetry generated by other mechanisms. On the other hand, thermal 
leptogenesis requires the reheat temperature after inflation, T R > M Rl pa 10 10 GeV. Such 
a high reheat temperature will generate a gravitino problem, unless the gravitinos are 
heavy with mass greater than O(10 TeV). 

Conclusions 

We construct a simple SO (10) SUSY GUT with D 3 family symmetry and an unbroken 
low energy R parity. The model describes fermion mass matrices with 14 parameters and 
gives excellent fits to 20 observable masses and mixing angles in both quark and lepton 
sectors, giving 6 predictions. Bi-large neutrino mixing is obtained with hierarchical 
quark and lepton Yukawa matrices; thus avoiding the possibility of large lepton flavor 
violation. The model naturally predicts small 1-3 neutrino mixing, sin $13 ~ 0.05, and a 
CP violating phase 5 close to tt/2. Among other interesting predictions is a tiny effective 
Majorana mass for neutrinoless double-beta decay. Leptogenesis is also possible with the 
decay of the lightest right-handed neutrino giving an acceptable CP violating asymmetry 
ei of order 10 -6 , and with the correct sign for the resultant baryon asymmetry. We also 
note that similar models with the non-abelian symmetry groups SU(2) or D^, instead 
of D 3 , can be constructed [T3|. 

The model presented in this letter provides an excellent benchmark for testing super- 
symmetric GUTs. At the very least, this model provides a phenomenological ansatz for 
fermion Yukawa matrices at the GUT scale, which fits low energy fermion masses and 

18 This formula is valid for standard model neutrinos. In the case of non-thermal leptogenesis in 
supersymmetric theories, the above value of ei is multiplied by an additional factor of 4, taking into 
account right-handed neutrino decays to leptons (or sleptons) in the final state with ordinary particles 
(and superpartners) in loops. 

19 In fact, we can obtain either sign. We obtain the opposite sign with a small increase of % 2 . 

20 Since we now also include the asymmetry due to sncutrino decays. 

21 S.R. would like to thank W. Buchmuller for emphasizing this point. 
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Table 1: The best fit for fermion masses and mixing angles. 

Initial parameters: (l/a G , M G , e 3 ) = (24.98, 2.77 x 10 16 GeV, -3.2 %), 

(A, e, a, e, p, e', f) = (0.64, 0.046, 0.83, 0.011, -0.053, -0.0036, 0.12), 

($ CT , $ e -, $ p , $ $ ) = (0.618, 0.411, 0.767, 3.673) rad, 

(m 16 , Mi/2, 4), A*(Af z )) = (3500, 450, -6888.3, 247.9) GeV, 

((m Hd /m 16 ) 2 , (m Hu /m 16 ) 2 , tan/3) = (2.00, 1.71, 49.98) 

(M^, M fi2 , M fll ) = (5.8 x 10 13 GeV, -9.3 x 10 11 GeV, 1.1 x 10 10 GeV) 



Observable 


Data ( <T ) 


Theory 


Pull 


Imfl^p^ in ripV l 

1 lllOiOOCO 111 VJ v~ V J 








n x i n 5 

w ^ /\ 1U 


1 1 fifiS7 CO 1 %) 

l.lUUiJI 1 U.l /C 1 


1 1fifi^5 


^0 01 

\ U.U1 




137 036 CO 1%) 


137.031 


< 0.01 


aJMr) 


1187 CO 002) 


0.1184 


0.02 


M 


174.3 (5.1) 


175.36 


0.04 


m b {M b ) 


4.25 (0.25) 


4.252 


< 0.01 


Mb — M c 


3.4 (0.2) 


3.513 


0.32 


m c (m c ) 


1.2 (0.2) 


1.03 


0.72 


m s 


0.105 (0.025) 


0.114 


0.13 


m d /m s 
Q- 2 x 10 3 


0.0521 (0.0067) 


0.0627 


2.53 


1.934 (0.334) 


1.763 


0.26 


M T 


1.777 (0.1%) 


1.777 


< 0.01 




0.10566 (0.1%) 


0.10566 


< 0.01 


M e x 10 3 


0.511 (0.1%) 


0.511 


< 0.01 


V us 


0.22 (0.0026) 


0.2192 


0.09 




0.0413 (0.0015) 


0.0407 


0.15 


v ub 


0.00367 (0.00047) 


0.00332 


0.56 


V td 


0.0082 (0.00082) 


0.00819 


< 0.01 


£k 


0.00228 (0.000228) 


0.00238 


0.19 


sin(2/5) 


0.736 (0.049) 


0.6757 


1.51 


Am§! x 10 3 


2.3 (0.6) 


2.407 


0.03 


Am^ x 10 5 


7.9 (0.6) 


7.866 


< 0.01 


sin 2 6 ) i2 


0.295 (0.045) 


0.2851 


0.05 


sin 2 6> 23 


0.51 (0.13) 


0.546 


0.08 


TOTAL x 2 






6.70 
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Table 2: Predictions for SUSY and Higgs spectra for the best fit given in Tabled 



Particle 


Mass (GeV) 


h 


121.6 


H 


495.5 


A 


499.9 


H+ 


411.5 


X? 


174.1 


xl 


241.7 


xt 


230.7 


9 


1181.0 


ti 


299.4 


h 


748.4 


fx 


1128.1 



Table 3: Predictions for neutrino masses, sin #13 and CP violation in the lepton sector 
for the best fit given in Table 



m V3 (eV) 


0.0492 


m U2 (eV) 


0.0096 


m vi (eV) 


0.0037 


sin 2 6»i3 


0.0024 


J 


0.0107 


sin 5 


0.98 


«i (rad) 


-1.286 


a 2 (rad) 


1.821 


(mpp) (eV) 


0.00021 


rrfJJ (eV) 


0.0065 


ei 


-1.16 x 10~ 7 
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mixing angles. There are several possible directions for future research which should be 
explored. First, the model will make predictions for flavor violating processes, such as 
\i — > e 7. This process, as well as other lepton flavor violating processes, will be tested 
to much higher accuracy in the future [IHj. Secondly, as discussed in the text we obtain 
acceptable values for the lepton number asymmetry parameter relevant for leptogenesis, 
ei . Thus it would be interesting to see if non-thermal leptogenesis via inflaton decay, as 
discussed by several authors (see |H] and references therein), can be incorporated into 
this theory. 

Finally, can this theory be derived from a more fundamental starting point such as 
string theory. Although we have not shown that this is possible, it is likely that it may 
be re- written as a five dimensional orbifold GUT (see [12] )• 
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